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We present our proposals for generating total hyperfine spin zero state for two / = 1 or two 
/ = 2 particles, starting from initial unentangled states. We show that our goal can be achieved 
by exploiting spin changing dynamics and quadratic Zeeman shifts with realistic choices of external 
magnetic fields and evolution time intervals. 
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I. INTRODUCTION 

State preparation for many-particle systems is gaining 
increasing attention recently. States with desired quan- 
tum entanglement properties find important applications 
in quantum information processing, cryptography, and 
computation Indeed efforts in this direction can be 
found in areas as diverse as quantum dots 0, trapped 
ions Q , and cold atoms 0, H| , to name a few. It has also 
been recognized recently that similar ideas can be used 
to create exotic many-body states in cold atoms. Many 
of these states are of particular interest in the condensed 
matter context. Examples of such proposals include the 
Haldane state for a spin-1 lattice @, resonant valence 
bond states Q, antiferromagnetic states Hf, d-wave 
"superconductivity" @, OH, Luttinger liquids [ll|, and 
even states that correspond to highest energy (instead of 
ground) state for a given Hamiltonian [6J, 1 1 21 ] . 

In this paper, we would like to discuss our proposal 
to generate a total hyperfine spin-zero singlet pair, start- 
ing from two particles that each has hyperfine spin-1 or 
2. Our primary motivation is that the ground states of 
certain one (spatial) dimensional lattice of Bosons with 
one atom per site and with suitable interaction between 
the Bosons are expected to be in the dimerized state 
[T3I [T3 ]. In this state, the lattice translational symme- 
try is spontaneously broken so that the spin-correlation 
between neighboring sites alternate between strong and 
weak. Furthermore, it can be shown that this dimer- 
ized state is smoothly connected to the ground state of 
a new Hamiltonian where alternate bonds are weakened 
[Til Il5j . In the limit where the alternate bonds are weak- 
ened completely to zero, we have a collection of isolated 
singlet pairs. The system remains gapped throughout 
this process. Reversing the argument, once singlet pairs 
can be prepared, in principle one can then obtain the 
many-body dimerized state of the original Hamiltonian 
by gradually turning on the coupling between the pairs. 
A similar argument can also be applied generally to va- 
lence bond states. Thus preparing the singlets would 
be a useful first step in generating these unconventional 
magnetic states. 

Irrespective of whether we can ultimately obtain the 
dimerized states experimentally using this method, we 
believe that our proposals of generating singlet pairs of 



hyperfine spin-1 and spin- 2 particles are of interest in 
their own right. While most literature on quantum infor- 
mation and cryptography deals with particles with only 
two internal states (e.g. spin-up versus down for spin- 
1/2 particles or horizontal versus vertical polarization for 
photons), it is known that higher dimensional quantum 
systems {i.e. systems with more than two internal states) 
offer many advantages such as higher tolerance in noise 
in quantum key distribution 16] or tests of Bell inequal- 
ities [ijj ■ A spin-singlet state of two spin-1 (2) particles 
is (one of) the most entangled two-particle state in a di- 
mension 3 (5) quantum system. It is therefore of high 
interest to know how to generate them, and here we of- 
fer a scheme starting from completely unentangled initial 
states 



We shall discuss two schemes. The first is simply 
to vary an external magnetic field adiabatically. This 
scheme is simple, but would work only in very special 
cases. The second, the central work of this paper, is 
a dynamical scheme. The main idea is to exploit the 
spin changing dynamics of two interacting particles in 
the presence of a (generally time-dependent) quadratic 
Zeeman field. This dynamics has been studied exper- 
imentally by Widera et al [19]. However, their focus 
was to extract the spin dependent interaction parame- 
ters (the spin dependent s-wave scattering lengths). We 
shall show that, for proper choice of parameters, one can 
prepare pure spin-singlet pairs using this spin changing 
dynamics. 



We introduce our proposal for hyperfine spin-1 parti- 
cles in section|TTJ The case for hyperfine spin-2 particles is 
discussed in section UTTl In these two sections, we empha- 
size the basic ideas and simplify our analysis by putting 
the magnetic field to be zero during some time intervals. 
The generalization to the cases where the magnetic fields 
are always finite are discussed in section HVl The conclu- 
sions are in section [V] In appendix A we provide some 
details of our estimate on the fidelity of our scheme. For 
simplicity, we shall occasionally use " spins" for hyperfine 
spins when there is no danger of confusion. 
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II. HYPERFINE SPIN-1 SYSTEMS 

We assume that one can initially prepare two particles 
in a potential well (or a double well, the detailed shape of 
the potential does not matter) , each with hyperfine state 
\f,m f ) = |/,0). In this (next) section / = 1 (/ = 2). 
This initial state has also been obtained by Widera et al 
(l9l [20| . We shall assume that these particles reside and 
stay in the lowest motional ground state throughout the 
dynamics, thus their spins are the only degrees of free- 
dom. A quantum mechanical state can therefore be indi- 
cated only by the quantum numbers m/ 's of the two par- 
ticles. For simplicity we shall suppress the / labels. Thus 
|0,0) m would indicate that both particles have mt — 0, 
whereas (|1, — l) m + | — 1, l) m ) indicates that there is 
one particle each in the m/ = ±1 state. Since we are 
dealing with Bosons, only symmetric states need to be 
kept. The subscripts m remind us that we are using the 
nif basis. 

Since we are interested in obtaining a total spin singlet, 
we would also like to keep track of the total hyperfine spin 
F of the two particles. It is therefore also convenient to 
use the basis \FMp) p where Mp is the projection of F. 
We only need to consider even F's. The subscripts F 
would be used to remind us when we are using this basis. 
The conservation of total m implies that we only need to 
consider the two states \00)f and |20)f- 

We consider the time evolution of these particles under 
the presence of a (in general time-dependent) quadratic 
Zeeman field. The linear Zeeman field needs not be con- 
sidered due to the conservation of total m and hence M 
within our proposed schemes. 

The Hamiltonian consists of simply two terms: the 
quadratic Zeeman energy and the interaction energy. 
The first contribution is a single particle term. It rep- 
resents am/ dependent energy quadratic in magnetic 
field B (for not too large _B's) and can be written in 

the form H Q = q fe mi , mj W + m%)\mi, m 2 ){mi, m 2 

where q — qB 2 with q a coefficient depending on the 
atom under consideration. The interaction energy is di- 
agonal in the total hyperfine spin basis \00)f and |20)f- 
We shall denote these interaction energies by Eq and E 2 
respectively. 

We first discuss the adiabatic scheme. The singlet state 
can be generated by adiabatically varying the external 
field B hence q if both the initial and final states are ei- 
ther the highest or the lowest energy states when q varies. 
Suppose q is positive. Then the initial state |0, 0) m is the 
lowest energy state for large q. If E2 > Eq, then the 
singlet state is also the ground state at q = 0. Hence, if 
the initial |0, 0) m state was produced at sufficiently large 
q, then decreasing the magnetic field to zero would auto- 
matically generate the singlet 1 00) p state. Such a scheme 
would not work of E 2 < Eq [2l[ , or their difference is too 
small for adiabaticity to be satisfied for q to be varied in 
time. 

Below we turn to our dynamical scheme which works 



for both E2 > (<)Eq. Let us first set up the Hamiltonian. 
In the |00) f and 1 20) f basis, the interaction part of the 
energy is simply a diagonal matrix, 



Hj 



E 
E 2 



Since we are targeting the state |00) f as the final state, 
we shall work in this basis. It is simple to rewrite Hq in 
this \FMp) basis since, from the Clebsch-Gordan coeffi- 
cients, we know that 



|00) F 
\20) F 



1 

1 



V6 

from which wc find 



(|l,-l) m + |-l,l) m -|0,0) m ) 
(|l,-l) m + |-l,l) m + 2|0,0) m ), (1) 



Q 3 



2 V2 
y/2 1 



Hence the total Hamiltonian is 



H = 



f E \ 




' 2 




{ E 2 J 


+!( 


v V2 


f) 



(2) 



In this basis, our initial state |0, 0) m is simply 
(w(0),v(0)) T = (-1/V3, v/273) T where T denotes the 
transpose. 

Let us consider the evolution of the state vector 
(u(t), v(t)) T in time t. Its value can be easily found since 
we have simply a two-level system. H can be rewritten as 
^ e +e 2 _|_ g\ +H e j f .f where r are the Pauli matrices and 
■A and HefLz - Eo^e* 



H, 



eff,x 



The scalar term 



only gives an overall phase and will be dropped below. 
The state vector at time t is then given by 



u{t) 
v(t) 



^cosfii 



H, 



eff 



-sinfii 



where 



n=\H eff \ = 



En — Ei 



9 



u(0) 
«(0) 



1/2 



,(3) 



(4) 



Here 2S7 is the Rabi frequency. The condition that we 
have a singlet at time t = t* is v(t*) = 0. Given that 
u(0) = -I/a/3 and v(0) = ^/2jl we need 



^psinm* 
and hence 



\ 2( costtf + 



cosfii* = 



sinf^ 



and 



H e ff,x + V2H e f 



= 0, 
(5) 

(6) 
(7) 
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These require, respectively, 



ir/fl 



and 



E 2 — Eq 



This last expression and Eq. (0} together give 



n = 



\E2 — Eq 



(8) 



(9) 



(10) 



Hence if q (hence q) has the same sign as E 2 ~ Eq, we 
obtain a singlet at times t* — |g ... As an example, 
consider 23 Na atoms in its / — 1 hyperfine state, where 
q > 0. The s-wave scattering lengths satisfy a 2 > a , 
hence E 2 > E . Hence it is possible to obtain a pure 
singlet at time t* when we choose the magnetic field to 

be5= ( E2 24 E " ) ■ The shortest possible time is t* = 
^2 = 2(&?-e ) " ^ spatial wavefunction is ^(r^), then 
E F = ^jjf fd 3 r\ip(r)\ 4 = a F U, where U is defined by 
this expression as in [l9j], and M a is the mass of the 
atoms. Immediately following that E 2 —Eq = U (a 2 — clq) , 
and for 23 Na a 2 — a is approximately 3.5a# }22j |. where 
is the Bohr radius. Given that q = 278 Hz/G 2 and if 
we choose Uas — (2-7r)30 Hz, then the required magnetic 
field and the time to obtain a singlet are calculated to be 
B = 1.089 G and t = 4.1 ms, respectively. 

As mentioned above this dynamical scheme is directly 
applicable to 23 Na, where q and E 2 — Eq are of the same 
sign. With simple modifications, this scheme can be gen- 
eralized to cases where q and E 2 — Eq are of the opposite 
sign, such as 87 Rb. This is done by simply letting the 
state first evolves under zero field till time t\, so that the 
wavefunction at t\ is 



1 p — iEgt-j 

V3 6 



2 p -iE 2 t 1 



e -iE 2 t! / _ e i<p 



V3 \ V2 



where <j> = (E 2 — Eo)ti is the phase difference between 
the F — and F = 2 states due to the interparticle 
interaction during the time interval < t < t\. Then a 
magnetic field and hence the quadratic Zeeman energy is 
turned on at ti, and the condition for obtaining a singlet 
after a further time t* (thus at total time t = t\ + t*) 
is simply the same as Eq. ([5]), except an extra factor e'^ 
multiplying the first term. In this case condition ([7]) is 
now replaced by 



H effyX cos<f) + V2H eff:Z = 0, 



which requires 



1 + 2cos(/> 



E 2 — Eq 



(11) 



(12) 



-(a) 



- 4 

3 
2 
1 



-(b) 





FIG. 1. Frames (a) and (b): two sets of solutions for ti and 
t* as functions of q. In frame (c), the value of -k/£1 is shown. 
Energy and time are in units of Eo — E2 and l/(Eo — E2), 
respectively. 



With suitable choice of cf) [i.e., ti), this condition can be 
satisfied even if q and E 2 — Eq are of opposite signs. For 
this q, the frequency £1 is given by 



n = 



\E 2 ~Eq\ 
1 + 2cos0| 



(2 



,1/2 



(13) 



and the required time t* is given by 



come 



' ' '^sinc 



V2n 

(sgng) 



(2 + cos 2 



a/2" 



(14) 



It follows that for a solution (cf)(ti),t*) solved by 
Eq. [|12H14|) , (0 + ?i27r,r +rmr/n) and (-0 + ?i2tt, -t* + 
mir/Q) are also solutions, where n and m are integers, 
and is an implicit function of (f>. Of the solutions above 
only positive t\ and t*'s are realistic, and Eq. (TT21 implies 
that there is no solution when — 1.5 < „ q P < 0.5 which 
holds for arbitrary signs of q and E 2 — Eq. Fig. [lja)(b) 
show two possible sets of (ti,t*), where E 2 < Eq and 
q > 0. In Fig. |Tfc) the value of ir/fl is shown. It is note 
that the solutions in Fig. [1] are also applicable for cases 
where q < and E 2 > Eq, due to the symmetry proper- 
ties of equations (TT2"|) and (TT4"|) . When the signs of q and 
E 2 — Eq are both reversed, it is observed that H — > —H, 
up to a constant, and <f> — 4> m Eq. (|T2|) . 

Consider 87 Rb atoms in its / = 1 hyperfine state, 
where the value of a 2 — ao is approximately — 1.4og p^| . 
By choosing Ua B = (2tt)30 Hz, then Eq - E 2 = 264Hz. 
If we further choose q = 1.5(Eq — E 2 ), then the required 
magnetic field is B = 2.35 G. Here we have used q = 71.7 
Hz/G 2 for 87 Rb. For this particular choice of q the solu- 
tions in Fig. HJa) and (b) are same and for the times are 
ti = 11.90 ms and t* = 3.435 ms. 

We shall estimate the error in producing the singlet 
state due to inaccuracies A^i and At* for t\ and t* , re- 
spectively. Let t\ and t* be the times needed to generate 
the singlet state so that \u(t)\ = 1 and \v(t)\ = 0, then the 
expected value of \u{M x ,At*)\ = y/l - KAii, Ai*)| 2 rj 

I 1 2 

1 — 4-. Here the approximation holds as long as Ati and 
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At* are not too large, and the analytical expression of 
\v(Ati, At*) | is shown in Appendix [A"l Take for instance, 
if | At! | + | At* | < 0.31 ms for 87 Rb, then \v\ 2 < 0.02 and 
|u| > 0.99. Similarly for the 23 Na, if the total inaccu- 
racies of time |Ati| + |Af*| < 0.36 ms, the same fidelity 
is achieved. We also note that from Appendix [3] and 
Fig- HJc), the value of \u(t)\ is closer to 1 for the same 
inaccuracies in time if a smaller value of \q\ is used. 

In the above, we have illustrated our scheme where 
zero magnetic field is utilized. In real experiments, zero 
magnetic field is not easily attainable, and undesirable 
effects may occur in low fields. To avoid those practi- 
cal issues, our scheme can be readily generalized to the 
case of finite field, by adjusting slightly the second ap- 
plied field, provided that the first field is not too large. 
Before we address this issue in section IIV1 we shall dis- 
cuss generalization of this dynamical scheme to / = 2 
systems. 



III. HYPERFINE SPIN-2 SYSTEMS 

The considerations here are direct generalization of 
those in the last section, and let us first consider the 
adiabatic scheme. Suppose q < 0, as in the case of / = 2 
hyperfine states of 23 Na and 87 Rb, then the initial state 
|0, > m is the highest energy level for large magnetic 
fields. In cases where Eq > i?2,4, then adiabatically 
switching off the external field would generate the de- 
sired singlet state. Conversely, if q > 0, the singlet state 
can be obtained if £2,4 > Eq. This scheme however is 
limited in its applicability, and in particular it neither 
works for 23 Na nor 87 Rb, where q < and 02,4 > ao- 

Let us proceed to consider the dynamical scheme. Sim- 
ilar to / = 1 spinors, we can write down the Hamiltonian 
for / = 2 spinors in the |00)f, |20)_f and |40)f basis, 
where 



|00>* = ^[|0,0>-(|1,-1> + |-1,1>) 
|20) F 



+ (|2,-2) + |-2,2))] T 
y-2|0,0) + (|l,-l) + |-1,1)) 

+ 2(|2,-2) + |-2,2))] 



(15) 



|40) F = ^[VT8|0, 0) + V8(\l, -1) + |-1, 1)) 



+ ^(|2,-2) + 1-2,2))] 



The Hamiltonian is 



H 



Eq 

E 2 
E 4 



2q 







/ 2 

14 17 12 
-12- 1 



,(16) 



where the two terms represent the interaction energy 
and the quadratic Zeeman effect, respectively. In 
this basis, our initial state |0,0) m is simply \ty ) = 
(l/-\/5, —y/2/7, ^/18/35) T , and the desired final state is 
the singlet state \00)f — (1,0, 0) T . In other words with 
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FIG. 2. In frame (a) the parameter regimes where we can not 
find solutions are shown as shaded regions. In frame (b), we 
show one solution (£1, t2, £3, £4) as function of q for 23 Na on 
the left and 87 Rb on the right. 



initial state \^o), we want to find a time-evolution oper- 
ator U so that 



F (20\U\* ) =f (40|L>|* ) = 0. 



(17) 



Here we need to solve for two complex equations, instead 
of just one (v(t) — 0) for / = 1 spinors. In the last 
section, we showed that our goal is achieved by evolving 
the system in two time intervals with a chosen q in each 
interval. Following this line, we attempt our goal with 
the following sequence: we first let the system evolve at 
q = for ti, and then at q ^ for t 2 , followed by q = 
for ta, then finally at q ^ again for Thus the time- 
evolution operator is 



U (ti, t2, t3, t4) 



e ~iHt 4 e -iH t 3 e -iHt 2 g — iffo*l 



u i u i u 2 u 1 



(18) 



Here Hq represents the Hamiltonian in Eq. (|16[) with no 
quadratic Zeeman effect (q = 0). 

It is convenient to define e(> 0) and 9 by the rela- 
tionships E2 — Eq = s cos 9 and E4 — Eq = esva.6. The 
amplitude e thus defines our energy scale and its inverse 
defines the time scale of the dynamics. The angle 9 then 
characterizes the element under consideration. For each 
(normalized) q and 9, we search for the times ti_4 so 
that Eq. (ITTl) is satisfied. We have limited ourselves to 
total time t\ + t 2 + £3 + <4 smaller than 200 in unit of 
1/e, so that the total experimental time scales remain 
reasonably short (see below) . There is a symmetry prop- 
erty: the solutions are the same for q = q* , 9 — 9* and 
q = —q*,9 = 9* + 7r (i.e., H — > —H, up to a constant). 

Fig. [U[a) shows the region where no solution is found 
(which may be due to the limits of our computational 
effort). Similar to the case of / = 1, when \q\ is small no 
solution is found. Here the values of 9 are limited to be- 
tween and 7r, given the symmetry property mentioned 
above. We can not find a solution at 9 = tt/2 (3tt/2) 
for q < (> 0). We note that for these special cases we 
allowed a longer search time with smaller stepsize, but 
no solution was found. When q approaches these values, 
fi_4 seem to diverge unfoundedly [231 ] . From the exper- 
imental point of view, if the value of 9 is very close to 
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these special cases, the system must have relatively long 
lifetime in order to achieve the singlet state using our 
scheme. However away from these values the singlet state 
can always be achieved with a reasonably short time. In 
Fig.[2ja), there also exist sharp spikes of q at 9 = 0, ir/4, 
tt/2, 7r and 3n/4 [24j. For general q's and 0's (excluding 
those discussed in the last paragraph), there exist many 
sets of solutions for £i_4, and no obvious relationship is 
found between different sets since the time evolution is 
not periodic (c.f. Sec HIj) . 

Let us consider 23 Na and 87 Rb in the / = 2 hypcrfmc 
state. For Na the differences of s-wave scatterin g le ngths 
(a>2 — o-o, a 4 — a>o) are approximately (11, 30)as [25j, and 
its quadratic Zeeman energy q = —278xB 2 Hz. For 87 Rb 
their values are (3.5, ll)as [25[ and —71.7 x B 2 Hz. Even 
though (d2 — ao, a 4 — ao) are different for 23 Na and 87 Rb, 
it happens that their #'s are similar, which are approx- 
imately 0.397T and 0.40-7T, respectively. As solutions U's 
arc not unique, wc only present one set of ii = i_ 4 as func- 
tions of g for both atoms in Fig.[^b). Take q/e = —0.5, 
then e • i l= i_ 4 =(2.645, 2.571, 5.555, 3.350) for 23 Na and 
(2.781, 2.442, 5.533, 3.343) for 87 Rb. Note that the dif- 
ference in 9 results in differences in t^s of up to 5%. In 
real experiments if we further choose Uas — (27r)30 Hz, 
then e = (2tt)959 Hz for Na and (2tt)346.5 Hz for 87 Rb. 
The controlled variables in laboratory units are B = 3.26 
G and t i= i_ 4 = (0.450, 0.437, 0.944, 0.570) ms for Na and 
3.90 G and (1.307, 1.148,2.601, 1.571) ms for 87 Rb. 

Let us also estimate the impact of inaccuracies in the 
time intervals used. With our estimation in Appendix [Al 
the expected amplitude in singlet state is more than 0.99 
if the total inaccuracies of time, |Aii| + |At2| + |Ai3| + 
|At 4 |, is less than 0.011 ms for 23 Na and 0.030 ms for 
8 Rb. Similar to / = 1 cases, to increase the tolerance 
in the time inaccuracies, we should use smaller value of 
\q\, with a caution that when \q\/e < 0.18 for both 23 Na 
and 87 Rb, no solution may be found. In other words the 
B must be larger than 1.95 G ( 23 Na) and 2.33 G ( 87 Rb) 
in order to arrive the singlet state. 

It is also possible to find solutions by restricting to 
three free time intervals, such as by letting t\ = 0, t 2 = 
t 4 , or t\ = t^, in the expense of fixing q to some special 
values. This is equivalent to the case in / = 1 systems, 
in which q = (E2 — Eq)/2 when we set t\ = 0. 



IV. NON-ZERO MAGNETIC FIELD RESULTS 

In the above sections we have shown the basic ideas of 
our proposal, exemplified with / = 1 and / = 2 spinors. 
However, zero B field environment may be difficult to 
prepare, and there could be undesirable effects in a zero 
field, such as spin-flips due to fluctuations of the magnetic 
field. Here we consider generalization of our scheme to 
the case where the magnetic field is always finite. For 
/ = 1, we let the system evolve at q = q\ 7^ for a time 
t\ and then at q = 02 7^ for a time t 2 . For / = 2, we 
first let the system evolve at q = q± =/= for a time interval 



ti, then at q = q 2 7^ for t 2 . followed by q = qi 7^ for 
another time £3, and finally again q = q 2 7^ for interval 
i 4 . For realistic solutions, q\ and q 2 must be of the same 
sign, dictated by the element under consideration. 



A. HYPERFINE SPIN-1 

Here we directly show the analytic solution for t\ and 
t 2 as functions of q\ and q 2 for / = 1. The state vector 
at time t = t\ + t 2 is given by 



u{t) 
v(t) 



cosQ 2 t 2 — i — e ^' 2 smil 2 t 2 



(19) 



x ( cos^i<i - i Hef ^ T sinfMi j ( j , 



where -ff e //,i and fii are defined as Eq. We obtain 
the solution of ti for the equation v(ti + t 2 ) = as 



<i = ±— cot 
Hi 



t 2 = =f— cot" 
\l 2 




T C(gi,g 2 ; 



C(gi,g 2 ; 



1/2N 



1/2N 



(20) 



where C(q[, q' 2 ) = 2q' 2 — '& — 6q' l — \2q' 1 q 2 . For convenience 
we have defined q- = qi/(E — E 2 ). From Eq. (|20| . given 
one solution (t\,t 2 ), the other solutions can be obtained 
as (t\ + nir/Vli,t 2 + to7t/£1 2 ) and (— t* + mr/f2i, — 1 2 + 
where n and m are integers. 
Interestingly no solution is found when the condition 
(2q' 1 + l)(2q' 2 + l){2q' 2 -3-6q' 1 -12q' 1 q 2 ) > does not hold. 
The conditions on q[ for the existence of solutions can be 
classified into two cases according to signs of q[; case(l) 
corresponds to q[,q' 2 < 0, and case(2) corresponds to 
q'i,Q 2 > 0. The conditions for both cases are listed below 
in more details. 



Case(l) 



Case(2) 



1 



-I<, 2 <!±^ 

2 - y2 - 2 - 12q[ 
or -^<ai<0, q' 2 < ~. 

1 3 3 + 6g; 

: < q < -, - < — < Qo. 

" qi ~ 6' 2 " 2- 12q[ ~ q2 



7<0 



Case(l) applies to 23 Na and case(2) applies to 87 Rb. 
When q[ = the results reduce to that in Section [H] 
In both cases either \q[\ or \q' 2 \ can not be too large, 
otherwise no solution is found. In Fig. |«3a) (b) we show 
the solutions for ti with the choices q[ — 1/96 and 1/12 
respectively for case (2), and solutions exist only when 
q' 2 > 49/30 and 7/2 accordingly. 

For the 87 Rb case, by choosing Uo,b = (27r)30 Hz, we 
have Eq — E 2 w 264 Hz. So the required magnetic fileds 
for q[ = 1/96 and 1/12 are B « 0.20 G and 0.55 G, 



6 




"(b) " 











2 4 6 „ c c x 4 6 4 6 

FIG. 3. In frame (a) and (b), plots ti and as function of q-2 
for two sets of solutions. In the frame (a) and (b) the values 
of qi/(Eo — E2) are 1/96 and 1/12 respectively. 
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FIG. 4. In frame(a) and (b) the parameter regimes where we 
can not find solutions are shown as shaded regions, where the 
values of qi/e are 0.002 and 0.02 respectively in the upper 
panels, and the values of qi/e are —0.002 and —0.02 respec- 
tively in the lower panels. 



respectively; for q' 2 = 49/30 and 7/2, B w 2.45 G and 
3.59 G, respectively. Solutions of (h,t 2 ) are (11.8,3.2) 
ms and (11.2,1.61) ms for (q'^q'z) = (1/96,49/30) and 
(1/12,7/2), respectively. 



the peaks near 9 = ir/4 and tt/2 grow quickly. Thus for 
< < 7r and q < (and similarly for ir < 9 < 2ir and 
q > 0), \q%\ cannot be too large in experiments. Besides 
there are also some changes in the peak positions, and 
some peaks are disappeared. 

For 23 Na and 87 Rb in / = 2, if q x /e = -0.002 and 
eti < 10, then solutions exist when q 2 /e —0.21. 
For example with (q 1 /e,q 2 /s) = (-0.002,-0.24) there 
is one solution of e • ti=i-4 = (3.76,4.56,3.72,2.38) for 
23 Na and (3.80,4.38,3.34,2.49) for 87 Rb. By choos- 
ing Uas = (27r)30 Hz, the B fields corresponding to 
{qi,q 2 ) are (0.20,2.28) G and (0.24,2.7) G, and the so- 
lutions for i i= i_ 4 are (0.624,0.757,0.618,0.395) ms and 
(1.75, 2.01, 1.54, 1.14) ms for Na and 87 Rb, respectively. 



V. CONCLUSION 

In this paper we propose a method to dynamically gen- 
erate a two-particle total hyperfine spin-singlet. For hy- 
pcrhnc spin-1 systems, by allowing the system to evolve 
in two time intervals, with a specific magnetic field in 
each interval. For hyperfine spin-2 systems, a similar 
scheme with four time intervals can also achieve our 
goal. For some special cases, such as negative q and 
E 2 ~ Eo < E±, the singlet state is however hard to ob- 
tain unless long evolution times are used. Except in these 
very special circumstances, our proposed scheme should 
be easy to implement in realistic experiments. Prepara- 
tion of these singlet pairs would allow us to create some 
exotic many-body states, and can also be useful in quan- 
tum information applications. 
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B. HYPERFINE SPIN-2 

For / = 2 spin systems in the non-zero magnetic field 
the time-evolution (|18|) operator is given as 

U{t 1 ,t 2 ,t 3 ,t 4 ) = e -iH 2 t Ae -iH lt3e -iH 2 t 2e -iH^^ ( 21 ) 

Here Hi represents the Hamiltonian in Eq. (|16l) . If q\ = 
0, the time-evolution operator is just Eq. (fT8|) . By a 
similar method in section IIIII the solutions ti 's can be 
obtained. Here we also have limited ourselves to Uot = 
^1+^2+^3+^4 < 200, in units of 1/e. There is a symmetry 
property: the solutions are the same when q^ = q* and 
6 = 6* or qi = -q\ and = -9*. 

In Fig. |U(a)(b) no solution is found in the shaded re- 
gions, where the values of qx/e are 0.002 and 0.02 for 
the upper-half of the panels, and the values of qi/e are 
—0.002 and —0.02 for the lower-half of the panels. Com- 
pared with the result of qi = (see Fig. HJ(a)) we find 
that when the absolute value of q± increases, the main 
change occurs in the region for negative q, where the 
total area with no solution found increases. Especially 
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Appendix A: Influence of inaccuracies in the time 
intervals used in evolution 

Here we provide some details on our error estimates on 
the final wavefunction. First we consider / = 1. With 
the condition v(ti,t 2 ) = 0, we can get the analytic form 
of \v{h + Ah,t 2 + At 2 )\ 2 ee | Aw| 2 as 

|Aw| 2 ~ ? p 20 Ah) 2 + (2g 2 At 2 f (Al) 

+ (AtiAt 2 )(12 ?l92 + 3E% - 6<?i£ 20 - 2^2^20)], 

where E 20 = E 2 — Eq. Here we have assumed that At\ 
and At 2 are small enough. In addition, Eq. (|A1[) can also 
be generalized to an inequality, 

|Aw| 2 < [fii|Ati| +n 2 \At 2 \f 

< [max(0 1 ,n 2 )(|Ati| + |Ai 2 |)] 2 = Al (A2) 
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To get this inequality the condition, (2q[ + l)(2q' 2 + 
l)(2q 2 — 3 — 6q[ — 12^2) > 0, has been used. In a 
special case where qi = 0, f^i reduces to |i?2o|/2 seen in 
section [TTJ 

This last formula for Ai gives an upper limit on the 
error produced due to the inaccuracies of time intervals 
employed in an experiment. If Ai = 0.02, then |u| (fi- 
delity) would be higher than 0.99. 

For / = 2 system we simply use the estimate |Aw| 2 < 



A 2 = [max(Awi,AiJ2)(|Ati| + |Af a | + |At 3 | + |At 4 |)] 2 , 
where Awi and Aw 2 are dehned as the maximum of 
\<*>i,i — w ij'l/2 and \u>2,i — W2j|/2, respectively. Here u>\ t i 
and u>2,i represent the eigenvalues of the Hamiltonian in 
Eq. (|16p at q = q\ and q = 52, respectively. Although 
we have not obtained the analytic form for the / = 2 
systems, we verified numerically that this formula for A2 
works well in practice. The fidelities discussed in text 
were estimated using this formula for A2. 
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